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Abstract 

We present a general method for constructing stochastic processes with prescribed 
local form. Such processes include variable amplitude multifractional Brownian mo- 
tion, multifractional a-stable processes, and multistable processes, that is processes 
that are locally a(t)-stable but where the stability index a{t) varies with t. In 
particular we construct multifractional multistable processes, where both the local 
<^ , self-similarity and stability indices vary. 

in 

^ : 1 Introduction 

O , In this paper we present a general framework for constructing stochastic processes with 
Q \ prescribed local forms. 

Stochastic processes where the local Holder regularity varies with a parameter t (usu- 
ally time) are important both in theory and in practical applications. The best known 
example is multifractional Brownian motion (mBm), where the Hurst index h of frac- 
tional Brownian motion is replaced by a functional parameter h{t), permitting the Holder 
exponent to vary in a prescribed manner. This allows local regularity and long range de- 
pendence to be decoupled to give sample paths that are both highly irregular and highly 
correlated, a useful feature for terrain or TCP traffic modeling. 

For modelling financial or medical data another feature is often important, namely the 
presence of jumps. Stable non-Gaussian processes give good models for data containing 
discontinuities, with the stability index a controlling the distribution of jumps. Recently, 
multifractional stable processes, generalising mBm, were introduced to provide jump pro- 
cesses with varying local regularity. However, a further step is needed for situations where 
both local regularity and jump intensity vary with time, for example to model financial 
data or epileptic episodes in EEG, where for some periods there may be only small jumps 
and at other instants very large ones. Our method may be used to construct processes 
where both h and a vary in a prescribed way: thus there are two parameters which might 
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correspond to distinct aspects of financial risk, to different sources of irregularity lead- 
ing to the onset of epilepsy, or to textured images where both Holder regularity and the 
distribution of discontinuities varies. 

It is natural to construct processes Y = {Y{t) : t G M} that have an identifiable local 
form near each u, that is where there is a limiting process 

y(.. + .r.)-yM ^ 

which may vary with u. If this limit exists as a non-trivial process we will say that Y is 
h-localisable at u and call the process Y^ = {Y^{t) : t G M} the local form of Y at u. The 
limit fll.ip may be taken in several ways: of particular interest are convergence in finite 
dimensional distributions, and convergence in distribtion; in the latter case we term the 
process strongly h-localisable. We will be especially concerned with /i-localisable processes 
with < h < 1 which are usually of a fractal nature. 

The most familiar example is mult if r actional Brownian motion Y which resembles 
index-h{u) fractional Brownian motion close to time u but where h{u) varies, that is 

, Y(u + rt) -Y(u) „ , , , , 

hm ^ L Li = 5 t 1.2 

where Bh is index-/i fractional Brownian motion, see [H [21 El [lOl [H]- Generalising this, 
mulitfractional a-stable processes have been constructed with local form /i(M)-self-similar 
linear a-stable motions |l9l EO] • 

It is clear from f ll.ip that the /i-local form at m, if it exists, must itself be /i-self- 
similar, that is Yl^{rt) = r^Y^(t) for r > 0. However, much more is true: under quite 
general conditions Y^ must be self-similar with stationary increments (sssi) at almost all 
u at which it is strongly localisable, that is r~^(Y^{u + rt) — Yliu)) = Y^(t) for all u 
and r > 0, see [H [9]. Thus if we wish to construct processes with given local forms, the 
local forms should themselves be sssi. Whilst this is a strong requirement, many classes 
of sssi processes are known, including fractional Brownian motion, linear fractional stable 
motion and a-stable Levy motion, see [6l ITT] . 

Our general construction will allow known localisable processes X{-,v) = {X{t,v) : 
t G M} for a range of v to be pieced together to yield a localisable 'diagonal' process 
Y = {X{t,t) : t G M} with local form depending on t. We will obtain conditions for 
the transference of the local properties of X{-,v) to Y. The basic setting is akin to that 
adopted in [21 US]. Thus we seek a random field {X{t,v) : G M^} such that for 

each V the local form X'^{-,v) of X{-,v) at v is the desired local form Y^ of Y at v. 
Typically, for each v the process {X{t,v) : t G M} will be one where the local form can 
be readily identified, such as an sssi process. Clearly the interplay of X{-,v) for f in a 
neighbourhood of u will be crucial to the local behaviour of Y near u. Thus the random 
field is set up as an integral or sum of functions that depend on t and v with respect to 
a single underlying random measure or process to provide the necessary correlations. In 
Section 4 we derive general criteria that guarantee the transference of localisability from 
the X{-,v) to y = {X(t,t) : t G M}; Section 5 addresses this for strong localisability. 

We illustrate the general method with several specific classes of processes. The method 
permits easy constructions of multifractional processes such as mult if r actional Brownian 
motion with variable amplitude (Section 6) and multifractional a-stable motions (Section 
7). In Section 9 we develop multistable processes, where the stability index a{t) is allowed 
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to vary. Here the constructions are based on sums over Poisson processes for which the 
required properties are reviewed in Section 8. In particular we construct multifractional 
multistable processes, where both the local self-similarity index and the stability index 
vary. 

2 Convergence of random processes 

This section is largely intended to establish notation. We work with two definitions of 
localisability of real valued random processes, one in terms of convergence of finite dimen- 
sional distributions and one requiring the stronger convergence in distribution, appropriate 
when the sample functions are viewed as members of some metric space. 

Given a probability space {^l,V, P), a random process X on a domain T is a family 
of random variables {X{t) : t G T}. For our purposes T will be either R or a subinterval 
of R, or sometimes a subset of R^ in which case we will refer to the process as a random 
field. 

We write — > X to mean that a family of random processes converges to a 
process X in finite-dimensional distributions. 

For processes with sample paths in suitable function spaces, convergence in distribution 
is defined in terms of a metric on the spaces. Let C (T) be the space of continuous functions 
on T C M. Writing d^{x,y) = sup^g^ — y{t)\ for the uniform metric on C{T), 

oo 

d{x,y) = Y,^-^^m{l,d^--'^\x,y)} {x,y e C{R)) (2.1) 

r=l 

defines a seperable metric on C(R) that gives the topology of uniform convergence on 
compact subsets of M. 

To accommodate processes with sample functions that have jumps, let T be a closed 
subinterval of R, and let D[T) denote the "cadlag" functions on T, that is functions 
which are continuous on the right and have left limits at all t E T. When T is a bounded 
closed interval we define a metric d^g'^^ on D[a, b] as follows. Let $ be the class of strictly 
increasing continuous bijections from [a,b] to itself. For each x,y E D[a,b] we define 
'^5 y) t)e the infimum of those 6 > for which there exists G $ such that both 
supo<t<i \4>{t) — t\<6 and supo<t<i \x{t) —y{(j){t))\ < 6. Then d^'''' is the Skorohod metric 
on D{[a,b]), see [15!, Chapter VI] or [3]. The Skorohod metric extends to a seperable 
metric on D(R) by 

oo 

ds{x,y) = J2^-^^m{l,d^-^'^\x,y)} {x,y e D{R)). (2.2) 

r=l 

Taking T as [a, b] or R, let F{T) be either C(T) or D{T) with the appropriate metric 
as above. Given a probability space {fl,V, P) we call X : Q ^ F{T) a random function 
or random element of F{T) if X^^{B) G V for every Borel subset B of the metric space 
F{T). If T' is a suitable subset of T and X is a random function on T then we may regard 
the restriction of X as a random function on T'. When we write X = F it will be clear 
from the context whether this refers to equality in finite dimensional distributions or in 
distribution. 
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For Xr and X random functions in F{T) where T is a closed interval, perhaps M, we 

say that Xj, converges in distribution to X, written X^ X, if E{f{Xr)) E(/(X)) 
for all bounded continuous / : F{T) M. Convergence in distribution is equivalent to 
convergence of finite dimensional distributions together with an appropriate stochastic 
equicontinuity condition, see |^ ITS]. 

Note that convergence in distribution in C(M) or D(]R) is equivalent to convergence 
in distribution of the restrictions of the random functions to every compact interval [a, b]. 
A technicality here is that our functions or processes may have a domain U that is a 
proper interval of R. This presents no difficulty, since Xr will generally be a sequence of 
enlargements of a process about some u interior to U, and the domain of definition will 
eventually include every bounded interval [a,b]. 



3 Localisable processes 



For convenience we give the definitions of localisability at u for random processes with 
domain R, but the definitions will also apply in the obvious way where the domain is 
a real interval with u as an interior point. Intuitively, a random process y on R is 
localisable at m G R if it has a unique non-trivial scaling limit at u. More precisely, we 
say that Y = {Y{t) : t G R} is h-localisable at u with local form the random process 
Y^a = {Y^it) : t G R}, if 

Y{u + rt) -Y{u) 



YL{t) (3.1) 

r" 

as r \ 0, where convergence is of finite dimensional distributions. If Y and have 
versions in C(R) or D(R) and convergence in (13.11) is in distribution, we say that Y 
is strongly localisable at u with strong local form Y^. Of course, strongly localisable 
processes are localisable, with the strong local form a version of the local form in C(R) or 
Z}(R). Note that the term locally asymptotically s elf- similar sometimes used for strong 
localisability. 

A number of well-known processes are /i-localisable, in particular processes that are 
h- self-similar, that is Y{rt) = r^Y{t) for r > 0, and which have stationary increments 
(that is Y{t -\- u) and Y{t) equal in law for m G R). 

Proposition 3.1 Let {Y{t) : t G R} be a process that is h- self- similar with stationary 
increments (h-sssi). Then Y is h-localisable at all u & with Y^ = Y . If in addition Y 
is in C(R) or D{R) then Y is strongly h-localisable at all m G R. 



Proof. If Y is /i-self-similar with stationary increments, then 
Y{u + rt) - Y{u) _ Y{rt) - Y{0) _ Y{rt) 



Yit) 



in law for all r 7^ 0, so F is localisible at u. 

Further, if Y is in C(R) or D{R) then {Y{u + rt) - Y{u))/r^ and Y{t) have identical 
probability distributions, since probability distributions on C(R) and D(R) are completely 
determined by their finite dimensional distributions, see [4]. Thus Y is strongly localisible. 
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There are several important processes which are sssi so which are strongly localisable 
by Proposition 3.1. 

For < h < 1, index- /i fractional Brownian motion (fBm) on M may be defined as a 
stochastic integral with respect to Wiener measure W: 

/oo 
f^t - xf-^l'' - {-xf;^!'') W{dx), (3.2) 
-oo 

where (a)+ = max{0, a} and ciK) is a normalising constant that ensures that the variance 
vari?/i(l) = 1. (Here, and throughout, we make the convention that expressions involving 
the difference of two positive parts represent an indicator function when the exponent is 
0, so for example, \ih = 1/2 then (t — a;)+^^^^ — (— a;)+^^^^ is taken to mean l[o,t)(x).) It 
is well-known [6], [7[ [131 E] that index-/;, fBm is an /i-self-similar process with a version 
in C(R) that has stationary increments, so is strongly localisable at all m G M with 

The a-stable processes form another important class of fractal processes of C(M), or 
of D{M) in the case of 'jump' processes, see Section 7. Under certain conditions a-stable 
processes may be sssi, see [171 Corollary 7.3.4], in which case by Proposition 3.1 they are 
strongly /i-localisable. 

A particular instance is linear stable fractional motion: 



J — oo 

(3.3) 

where < a < 2 and M is an a-stable random measure with constant skewness (3 and 
control measure Lebesgue measure, < h < 1 and a and b are constants, see [T71 Section 
7.4 and Chapter 10]. The process is /i-sssi and so is /i-localisable at all m G M with 
{La^hYu = La^h- Provided that h > 1/a it has a version in C(M), so is strongly localisable. 
However, if h < 1/a then almost surely Y is unbounded on every interval and so is not a 
process of -D(ffi), though it is nevertheless localisable. (Note that later we will represent 
such processes as Poisson sums rather than integrals with respect to random measures.) 

An a-stable Levy motion, < a < 2 is a process in D{M.) with stationary independent 
increments which have a strictly a-stable distribution. It may be represented as 

L„(t) = M([0,t]) (3.4) 

where M is an a-stable random measure on M with constant skewness intensity, see [T71 
Section 7.5]. Then La is 1/a-sssi, and so is strongly l/a-localisable. 

In later sections we will give general constructions of localisable processes where the 
local form varies with u. For now we note that localisability behaves well under 
reasonably smooth changes of coordinates. In particular the following proposition allows 
the introduction of varying 'local amplitude' for localisable processes. 

Proposition 3.2 Let U he an interval with u an interior point. Suppose that {Y{t) : t G 
U} is h-localisahle (resp. strongly h-localisahle) at u. Let a : U ^M. satisfy an rj-Holder 
condition on U, that is 

|a(t) -a(t')| < clt-t'l" {t,t'eU), 

where rj > h. Then aY = {a{t)Y{t) : t G f/} is h-localisahle (resp. strongly h-localisahle) 
at u with {oY)'^ = a{u)Y^. 
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Proof. We have 

a(u + rt)Y (u + rt) — a(u)Y (u) , , F(m + rt) — , afw + rt) — afw) 
^ — T-^ \ I \ I ^ _^ J w _^ ^ L^. 

The result now follows on letting r — > with the appropriate form of convergence, noting 
that the right-hand term has zero limit almost surely. ■ 



4 Localisable processes with prescribed local form 

We aim to construct localisable funtions with prescribed local form by 'joining together' 
locahsable processes {X(t,v) : t e U} over a range of v. Thus we seek conditions that 
ensure Y — {X(t, t) : t E U} looks locally like {X{t, u) : t E U} when t is close to u. 

Let U be an interval with u an interior point. Let {X{t,v) : {t,v) E U x U} be a 
random field and let Y be the diagonal process Y = {X{t,t) : t E U}. We want Y and 
X{-,u) to have the same local forms at u, that is Y^{-) = X!^{-,u) where X'^{-,u) is the 
local form of X(-, m) at u. Thus we require 

X{u + rt,u + rt) - X{u,u) idd s 

-h >Kit,u) (4.1) 

as r \ 0. The following theorem gives a sufficient condition for this to occur. 

Theorem 4.1 Let U be an interval with u an interior point. Suppose that for some 
< h < r) the process {X{t, u) : t E U} is h-localisable at u E U with local form X!^{-, u) 
and 

P{\X{v,v)- X{v,u)\ > Iv-u]"^) (4.2) 

as V ^ u. Then Y = {X(t, t) : t E U} is h-localisable at u with Y^{-) — X^{-, u). 
In particular, this conclusion holds if for some p > and rj > h 

E{\X{v, v) - X{v, u)\P) = 0(1^; - ul"^) (4.3) 

as V —>■ u. 

Proof. For r 7^ 

Y{u + rt)-Y{u) _ X{u + rt,u + rt) - X{u,u) 

X{u-\-rt,u-\-rt) — X{u-\-rt,u) X{u + rt,u) — X{u,u) 

Fix t eM. and c > 0. Let ro be sufficiently small to ensure that if < r < ro then both 
u±rt eU and cr'* > {r\t\)''. Then for < r < Tq 

p , \X{u + rt,u + rt) -X{u + rt, u) \ ^ 

< P{\X{u + rt,u + rt) -X{u + rt,u)\ > (r|t|)'?) 

< P{\X{u + rt,u + rt) - X{u + rt,u)\ > \{u + rt) - u]"^) 
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as r \ 0, by Thus for all t e R, 



X{u + rt,u + rt) — X{u + rt, u) 




since is localisable at u. We conclude from (14. 4p that Y is localisable at u with 

local form FJ(-) =X'^{-,u). 

If (14. 3 p holds, Markov's inequality implies (14. 2 p (with rj replaced by some h < rj' < t]) 
and the conclusion follows. ■ 

Although Theorem 4.1 is valid for all /i > 0, it is normally applied with < h < 1. 
If X{-,u) is /i-localisable for h > 1 then the limit of (14.40 is usually dominated by the 
left-hand term giving that Y is 1-localisable, see Theorem 9.4 for an example of this. 

5 Strongly localisable processes with prescribed local 
form 

We obtain an analogue of Theorem 4.1 in the strongly localisable case, that is a criterion 
for convergence in distribution in (13.11) . 

Theorem 5.1 Let F(R) be either C(M) endowed with the metric d or D(M) with ds, see 
(\2.1\ii or ( \2.i^ . Let U he an interval with u an interior point. Suppose that for some h > 
the process {X{t, u) : t E U} of F{U) is strongly h-localisable at u, with local form X4(-, u) 
a random function of F{M.). Suppose that for all c> 



as e ^ 0. If the process Y = {X(t,t) : t E U} is in F{U) then Y is strongly h-localisable 
at u with Y^{-) = X'^{-,u). 

In particular, this conclusion holds if for some t] > h we have 




(5.1) 



sup 



X{y,v)- X{v,u)\ 



< oo 



(5.2) 



V — uy^ 



almost surely. 



Proof. First consider (C(]R),(i). For each positive r and r sufficiently small. 




(5.3) 
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X(u + rt,u + rt) — X(u + rt,u) 
as r — s> 0. Thus, the restriction oi to [— r, rj, converges to 

in probabiUty in (C[— r, r], d'"'^''^') as r — > 0. From the definition (12.1 p of d, convergence 
in probabihty on every bounded interval imphes convergence in probabihty on (C(]R), d), 
so 

X{u + rt,u + rt) — X{u + rt, u) p 



(5.4) 



in {C{M),d). Then 

Y{u + rt)-Y{u) X{u + rt,u + rt) - X{u,u) 



X{u + rt, M + rt) — X{u + rt, u) X{u + rt, m) — X(u, m 



+ 



^ X'uit,u), (5.5) 

as r — i> 0, since X is locahsable at m. Here we use a standard property [H Theorem 
4.1], that for random elements Z,., Z, Wr, W of some metric space {M,p), ii Z^. Z and 

p(Z,, W,) ^ 0, then ly, ^ Z. 

Turning to {D{^),ds), if X(t, n) G D(]R) and (15. ip holds, the same argument using 
(15. 3p implies convergence in probability in (15. 4p with respect to the metric ds- (Note that 
(^g'^''^\f, 0) < supjgj_^,_] |/(t)| for each r for / G -D(R).) Convergence in distribution then 
follows just as for C*(M). 

Finally, (15. ip is an immediate consequence of (15. 2p if h < rj. m 

To utilise Theorem 5.1 we need to verify (15. 2p . that is to show that Z{v) = {X{v, v) — 
X{v,u))/\v — is bounded as v ranges across an interval. The following form of Kol- 
mogorov's continuity theorem will be extremely useful for this. 



Theorem 5.2 (Kolmogorov's continuity theorem) Let {Z{v) : v G T} be a random pro- 
cess where T is a bounded subset o/M". If for some p > 0,e > and c > 

E\Z{v) - Z[v')Y < c\v - {v,v' G T), 

then Z has a continuous version that is almost surely rj-Holder continuous for all < rj < 
e/p. 



Proof. See, for example, [HI Theorem 25.2]. ■ 



6 Multifractional Brownian motion with variable am- 
plitude 

A number of constructions of multifractional Brownian motion, a process with index- 
h{u) fractional Brownian motion as its local form at u, have been given, see [H [21 El 
To demonstrate our method we indicate briefiy a straightforward construction of 
multifractional Brownian motion, that is strongly localisable with a given local index and 
amplitude. 
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As in p!^ we model our definition on (13.21) but allow h to vary. By virtue of Proposition 
3.2 variable local amplitude presents no difficulty. Let U he a bounded closed interval 
and let /i : — > (0, 1) satisfy an 77- Holder condition 

\h{v)-h{v')\<k\v-v'\'^ {v,v'eU) (6.1) 

where < rj < 1. 

Theorem 6.1 (Multifractional Brownian motion) Letu G M and let U he a closed interval 
with u an interior point. Suppose that h : U ^ (0, 1) and a : U ^ both satisfy an 
Tj-Holder condition where h{u) < rj <1. Define 

/oo 
({t - x)f - {-x)f^-^n W{dx) {t G U). (6.2) 
00 

Then Y is strongly h{u)-localisahle at u with = a{u)c{h{u))Bh(^u) where is index-h 
fBm and where c{h) is the normalisation constant in { \3.^ . 

Proof. By Proposition 3.2 it is enough to consider the case where a{v) = 1. We define a 
random field by the stochastic integral 

/oo 
Ut-x)'i''^-''^ -{-x)f^-''^)w{dx) {t,veU), (6.3) 
-00 ^ 

where W is Wiener measure on M. Since the integrand of (16. 3p is square integrable, X{t, v) 
exists a.s. with mean for all t,v & U. A mean value estimate applied to (16.31) easily 
gives that E(X(t, v) - X{t', v')y < c{\t - t'l^"' + \v - y'l^"^), first for p = 2 and then for all 
p > since the increments are Gaussian. By Kolmogorov's criterion, for all e > such 
that h{u) < rj — e there is an a.s. finite random variable C such that 

\X{v,v)-X{v,u)\<C\v-u\'^-' (veU) 

so (15. 2p holds (with rj replaced by r/ — e). But X{-,u) = Bh[u){-) which is sssi so is strongly 
/i(M)-localisable at u by Proposition 3.1. Theorem 5.1 implies that Y = {X{t,t) : t G T} 
is strongly /;,(M)-localisable at u with yj(-) = X'^{-,u) = {Bh(.^))'^{-) = Bh(u){-)- ■ 

7 Multifractional stable processes 

Multifractional Brownian motion generalizes fractional Brownian motion by allowing the 
parameter h to vary with time. By working with a stochastic integral with respect to 
an a-stable measure instead of Wiener measure, we now construct multifractional stable 
processes with the local scaling exponent depending on t. 

Recall that a process {X{t) : t G T}, where T is generally a subinterval of M, is 
called a-stable (0 < a < 2) if all its finite-dimensional distributions are a-stable, see 
the encyclopaedic work on stable processes [Uj. Note that 2-stable processes are just 
Gaussian processes. 

Many stable processes admit a stochastic integral representation. Write Sa{cr,P,fi) 
for the a-stable distribution with scale parameter a, skewness /? and shift-parameter /i; 
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we will assume throughout that /i = 0. Let {E,S,m) be a sigma-finite measure space 
(which will be Lebesgue measure in our examples). Taking m as the control measure and 
(3 : E ^ [—1,1] a measurable function, this defines an a-stable random measure M on E 
such that for A e ^ we have that M{A) ~ ^„(m(A)i/", J^/3{x)m{dx)/m{A), 0). If /? = 
then the process is symmetric ct-stable or SaS. 
Let 

J^a = ^aiE,£,m) = {/ : / is measurable and < oo}, 
where || \\a is the quasinorm (or norm if 1 < a < 2) given by 

{J^\fixrmidx)f'' (a^l) ^7^^) 

Ie \f{x)\m{dx) + \f {x)P{x) In \f{x)\\m{dx) {a = 1) 

The stochastic integral of / G J^aiE,S,m) with respect to M then exists [T71 Chapter 3] 
with 

Hf) = f f{x)M{dx) ~ S^{af,[3f,0), (7.2) 
Je 

where 

Jf{x)<'^>P{x)m{dx) 

-/-ll/IU Pf- , 

writing a^'^^ = sign(a)|a|'', see [T71 Section 3.4]. In particular, 

E|JC/■^|P = / c(a'/^'P)ll/lla {0<p<a) ,^ 
\ )\ ^ oo (p > a) ^ ' ^ 

where c{a,P,p) < oo, see [171 Property L2.17]. 

When < a < 1 there is a non-negative stable subordinator measure M' associated 
with M so that M'{A) ~ S'«(m(A)i/°, 1, 0). In particular, for / e J^a, 



|/(/)| < / \fix)\M'{dx). (7.4) 

We will be concerned with processes that may be expressed as stochastic integrals 

X{t)= I f{t,x)M{dx) + fx{t), (teT), (7.5) 
Je 

where f{t, ■) is a jointly measurable family of functions in J^a{E,S,m) and /i(t) are real 
numbers. Note that if esssup^<4<^/(t, x) = oo for all x G A for some A C E with m{A) > 
then X{t) will be unbounded a.s. on the interval [a,b], see [T71 Section 10]. 

Here we consider the localisibility at u of processes defined in terms of random fields 

X{t,v)= / f{t,v,x)M{dx) + fi{t,v) {t,veU) (7.6) 
Je 

where f(t,v,.) G J-'a and fi(t,v) G M for all t,v & U for some interval U. We assume 
throughout that f(t,v,x) is measurable onU x U x E. 

The term /i(t, f ) is easily dealt with: if v ^ fi{v, v) is pointwise r^-Holder at f = m, 
that is \fi{v,v) — fi{u,u)\ < k\u — for v close to u, where < h < rj < 1, then the 
/i-localisability of y = {X{t,t) : t E U} at u and its local form are unaffected if we set 
li{t,v) = 0, so we assume this throughout this section. 

The following proposition gives conditions for Y to have a continuous or bounded ver- 
sion, which is needed for strong localisability to be meaningful. Note that these sufficient 
conditions are geared towards our context; for other aspects see [TTl Chapters 10,12]. 
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Proposition 7.1 Let U be a closed interval. Let X be a random field defined by 

X{t,v)= [ f{t,v,x)M{dx) {t,veU) (7.7) 



where f{t,v, ■) G J-'a are jointly measurable and M is an a-stable random measure with 
control measure m and measurable skewness. 
(a) Let < a < 1. // 

sup \ f{t,v,x)\\\a < oo, (7.8) 



then the random field ( [7. 7[ ) has a bounded version. 

If in addition {f{t,v,x) : x G E} is an equiuniformly continuous family for t,v & U, 
then (|7. ?P has a continuous version. 



(b) Let 1 < a <2 and 1/a < r] < 1. If 

||/(t, V, ■) - fit', v', .)\\^<k{\v-vr+\t- tT) (t, t', V, v'eU), (7.9) 

then Y = {X{t,t) : t G f/} has a continuous version for t G U, satisfying an a.s. (3 -Holder 
condition for all < f3 < {rja — l)/a. 

Proof, (a) Since < a < 1 there exists a stable subordinator measure M' associated 
with M, so that M' has control measure m and M'{A) ~ Sa{m{AY/°', 1, 0). By (El), for 
t,v eU, 



\X{t,v)\< / \f{t,v,x)\M'{dx) < / sup \f{t,v,x)\M'{dx) = Z, 
J J t,veu 

where Z is an almost surely finite random variable by (17.81) . so X(t, v) is a.s. bounded for 

t,v eU. 

Now assume also the equicontinuity condition. Given e > we may, since E is a- 
finite, choose D C E such that /g^^ (sup^^^g^ \f{t,v,x)\)°'m{dx) < e". By equiuniform 
continuity we may find 6 > such that for all x & E and \{t,v) — {t',v')\ < S we have 
\fit,v,x) - f{t',v',x)\ < m(L>)-i/"e. Then if \{t,v) - {t\v')\ < 5, 1^ gives 

\X{t,v)-Xit',v')\ < [ \f{t,v,x)-f{t',v',x)\M'{dx) 



E 



< 2 



E\D t,veU 



sup \f{t,v,x)\M'{dx)+ [ ' M\dx)^Z, 



say, where Z^ is a random variable. Fix < p < a. By (17. 3p there is a constant c 
independent of e such that 

E\Z,\P< ceP. 

Thus choosing e(n) {n = 1,2, . . .) such that E|Ze(„)p < 2~", there are corresponding 6n 
such that 

sup \X{t,v)-X{t',v')\<Z,^r.). 

\it,v)-{t' ,v')\<Sn 

Since Yl'^=i E|^e{n)|'' < oo, the Borel-Cantelli lemma gives that Z^(^n) almost surely, 
so sup|(t \X{t,v) — X{t',v')\ — >• a.s. as n — >■ oo, giving continuity of X{t,v) 
a.s. 
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(b) From 

X{t,v)-X{t',v') = j {f{t,v,x)-f{t',v',x))M{dx). 

This integrand is in JF^,, so for < p < a, estimate (17.31) gives 

E\X{t,v)-X{t',v'W < c^\\f{t,v,-)-f{t',v',-)t 

by (17.91) where Ci and C2 are independent of t, t', v, v' E U. Speciahsing, 
E\Y{t) -Y{t')\P = E\X{t,t) - X{t',t')\P < 2c2\t-t'\''P 

for t, t' G U. 

Since t] > 1/a we may choose < J9 < a such that rjp > 1. Kolmogorov's Theorem 
5.2 gives that Y has a continuous version for t E U with an a.s. /3-Holder condition for 
all < < {rjp — l)/p for all p < a. m 

We require the following calculus lemma. 

Lemma 7.2 Let U he an interval and let f : U ^ W be continuously differentiate with 
f satisfying an rj-Holder condition 

\f'{v)- f'{w)\<k\v-w\'^ {v,weU) (7.10) 

for some < rj < 1. Let v,w,ueU with v ^ u,w u. Then 

f{v)-f{u) f{w)-f{u) 



V — u w — u 



< 2U|t;-u;|''. (7.11) 



Proof. We may assume without loss of generality that v < w and u < w. Write 

y/ \ f 

g(v) = . We consider three cases. 

V — u 

(a) If f < M < w, then by the mean value theorem there exist vq G {v,u) and 
Wo G {u,w) such that g{v) = f'{vo) and g{w) = f'{wo). Then 

\giv) - g{w)\ = \f\vo) - fiwo)] < k\vo - wo^ < k\v - w\^. 

(h) li u < V < w and |w — f| > |f — m|, then — f| > lliu — m|. There exist vq G (u, v) 
and Wo G {u,w) such that g{v) = f'iyo) and g{w) = f'{wo), so 

\g{v) - g{w)\ = \f'{vo) - /'(«;o)| < k\vo - wo\'' < k\w - u\'' < kr'\w - v\^. 

(c) li u < V < w and |w — f | < |f — u|, we apply the mean value theorem to g. Thus 
there exists s G {v,w) such that 

g{v)-g{w) = {v-w)g'{s) 

(s-u)f'is)-f{s) + fiu) 



{v — w) 
(f — w) 



(s — uy 

f'is) - f'iz) 



[s - u 
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\v — 


w\ 


\s - 


- z\ 






\s - 


- u\ 





where z e {u, s) using the mean value theorem again. Hence 

\g{v)-g{w)\ < k 

< k\v — wWs — u\'^~^ 

< k\v — w\\v — ul'^"^ 

< k\v — w\'^. 



The following theorem gives conditions that allow the transfer of localisability prop- 
erties from X{-,u) to F = {X{t,t) : t G U} in the a-stable case, generalising the results 
of Section 6 in the Gaussian case. 

Theorem 7.3 Let U be a closed interval with u an interior point. Let X be a random 

field defined by 

X{t,v) = J f{t,v,x)M{dx) {t,veU) (7.12) 

where f{t,v,-) e o-f^ jointly measurable and M is an a-stable random measure with 
control measure m and measurable skewness. 

(a) Suppose that < a < 2 and the process X{-,u) is h-localisable at u with h > 0. 
Suppose that for some r] > h 

\\f{t,v,-)-fit,u,-)\\a<ki\v-u\^ {t,veU). (7.13) 

Then Y = {X{t,t) : t E U} is h-localisable at u with local form Y^^-) = X'^{-,u). 

(6) Suppose that < a < 1 and that X{-,u) is strongly h-localisable in C(R) (resp. 
D{M.) ) at u. Suppose that for some rj > h 

\f{t,v,x)- f{t,u,x)\<ki{x)\v-u\'' {t,veU,xeE), (7.14) 

where ki{-) G J-'a- If Y = {X(t,t) : t G U} has a version in C{U) (resp. D{U)) (see 
Proposition 7.1(a)), then Y is strongly h-localisable at u in C(M) (resp. D{R)) with 
Yi{.)=X'S-,u). 

(c) Suppose that 1 < a < 2, that rj > 1/a and that X{-,u) is strongly h-localisable in 
C(]R) or I?(M) at u. Suppose that for all t,v & U the partial derivative fv{t,v, •) G 
with 

\f,{t,v,x)- f,{t,v',x)\<ki{t,x)\v-v'\'> {t,v,v' eU.xeE), (7.15) 
where sup^^j; ||A;i(t, ■)\\a < oo, and that 

sup\f^{t,v,x)-f^{t',v,x)\<k2{t,t',x) {t,t' eU,xeE), (7.16) 

where \\k2{t,t\ < c\t — t'\^K Then Y — {X{t,t) : t & U} is strongly h-localisable at u 
mC{^) wtthYi{-) = X'^{;u). 

Proof, (a) We have 

X{t, v) - X{t, u) = j {fit, V, x) - fit, u, x)) Midx) (7.17) 
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so, taking < p < a and using (17.31) . there is a constant ci such that 

E\Xit,v)-X{t,u)\^ < c,\\fit,v,-)-fit,u,-)\\l 

< cih\v-u\''P. 

The conclusion follows from Theorem 4.1. 

(b) Since < a < 1 there exists a stable subordinator measure M' associated with 
M, so that M' has control measure m and M'{A) ~ ^^(m(A)i/", 1,0). Applying to 
fl7.17p and using (17.141) . gives that for t,v e U 

\X{t,v)-X{t,u)\ < [\f{t,v,x)-f{t,u,x)\M'{dx) 



< \v-u\'' j ki{x)M\dx) 

< \v-u[^Z, 

where Z is an a.s. finite random variable. Thus (15. 2p holds and Theorem 5.1 gives that 
Y is strongly localisable at u. 

(c) It is easy to check that Y satisfies the conditions of Theorem 7.1(b) and so has a 
continuous version. Again we verify (15.21) . Define 

Z(t,v) = ^ ' ' (t,veU,Vy^u) 

V — u 

g(t, V, x)M{dx) 
where 

f{t,v,x) - f{t,u,x) 

g{t,v,x) = . 

V — u 

Applying Lemma 7.2 with /(t>) = f{t,v,x) and noting (I7.15p . we get 

\g{t,v,x) - g{t,v,x)\ < 2'^ki{t,x)\v - v'l"^. (7.18) 

Also 

\g{t,v,x) -g{t',v,x)\ = ^^3^ \{f{t,v,x) - f{t',v,x)) - {f{t,u,x) - f{t',u,x))\ 

< \fy{t,vux) - f^{t',vi,x)\ 

< k2{t,t',x), (7.19) 

for some Vi G {u,v), on applying the mean value theorem to f{t,v,x) — f{t',v,x). From 
(I7.18P and (17.190 together with the conditions on ki and k2 we get 

\\g{t, V, ■) - git', v', ■)L<c,i\v- vT + \t - t?). 

Applying Proposition 7.1(b) to {Z{v,v) : v G U}, it follows that Z{v,v) = {X{v,v) — 
X{v, u))/{v — u) has a version that is a.s. continuous and bounded for v eU . Thus (15. 2p 
holds and strong localisability follows from Theorem 5.1. ■ 
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We illustrate Theorem 7.3 by constructing processes whose local forms are linear stable 
fractional motions La,h{t), see (13. 3p . Overlapping results with a different emphasis are 
given in [191 [20]. The following process is termed a linear stable multifractional motion: 

a ((t - x)f - (-x)f )-^/") 

oo 

+ 5 (^(t _ a;)^W-V" _ j j (^teR), (7.20) 

where M is an a-stable random measure (0 < a < 2) with constant skewness intensity (3 
and control measure Lebesgue measure, with h(t) G (0, 1) for all t G M, and a and b real 
numbers. (Recall that (w)+ = max{0,?i;} and (w)^ = — (w)+ for w G M.) 
To investigate localisability, we introduce the random field 

a((t-x)^")-'/"-(-x)f)-'/" 

-oo 

+ 6 ((t - x)'!^''^-'/" - (-x)^('^)"'/") ] M(tix) {t,veR). (7.21) 

Then X{t, v) is well-defined since since for each {t, v) the a-th power of the integrand is 
Lebesgue integrable. For each fixed v the process X{-,v) is just a linear stable fractional 
motion ( 13. 3p so is /i(f )-localisable, with X'^{-,v) = La,h(v){-) for all m G M. Provided that 
h{v) > 1 /a it is in C(R) and is strongly localisable. 

Theorem 7.4 (Linear multifractional stable motion) Let U be a closed interval with u 
an interior point. Let < a < 2 and h : U — > (0, 1). Define {Y{t) : t e U} by (\7.20\j . 

(a) Assume that h satisfies a rj-Holder condition at u 

\h{y) - h{u)\ < k\v - {v e U) 

where h{u) <rj <1. Then Y is h{u) -localisable at u with local form Y^ = L^^^u)- 

(b) If 1 < a < 2 and h is differentiable with 1/a < h{u) < 1 and 

\h\v)-h\v')\<k\v-v'\^ {v,v'eU) (7.22) 

where 1/a < rj <1, then Y is strongly h{u) -localisable at u with local form Y^ = L^t^hiu)- 

Proof. For brevity of exposition we give the proof in the case of well-balanced linear 
multifractional stable motion, that is with a = 6 = 1 in (I7.20p and (I7.2ip : the general 
case is very similar. Thus we take 

f{t,v,x) = |t - a;|^(^)'i/" - |a;|^(^)-i/" 

in Theorem 7.3 (when h{v) = 1/a such expressions are interpreted as l[o,i](x) where l[o,t] 
is an indicator function). Then X{t,v) = J f{t,v,x)M{dx) and Y{t) = J f{t,t,x)M{dx). 

(a) By continuity, we may assume that [/ is a sufficiently small interval to ensure that 
h{v) < f] for all f G f/. Fix h_, such that < /i_ < h{v) < < 1 for all v eU . Then 
for each t, f , v',x E U with x ^ 0,x ^ t, the mean value theorem gives 

\f{t,v,x) - f{t,u,x)\ 

= ||t-x|^(-)-i/°log|t-x| - Ixl'^^-^-^/^loglxll \h{v)-h{u)\ 

< ||t-a;|^(-)-i/°log|t-x| - |x|'^(-)-i/"log|x|| A;|f (7.23) 
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where h{-) = h{t,v,x) G [h{v),h{u)]. But 

k\\t- xl'^^-)-^/" log \t-x\- Ixl'^^-)-^/" log |x| I < ki{t, x) 
for all t E U,x E M, where 



cimaxjl, + < 1 + 2 ma^t^u \t\) 

(|x| > l + 2maxt6t/|t|) 



for appropriately chosen constants ci and C2- Then J ki{t,x)°'dx is finite and uniformly 
bounded for t E U, so as is /i(M)-locahsable at u, Theorem 7.3(a) gives that 

Y = {X(t,t) : t E U} is /i(M)-localisable at u with local form Y^{-) = X^(-,m) = La,h(u){-)- 
(b) We may assume that U is small enough and h-, /i+ are chosen so that < 1/a < 
h_ < h{v) < < 1 for all v E U. A similar estimate to fl7.23p on the derivatives gives 

\Mt,v,x)-Ut,v',x)\ < [\h'iv)\\h\vi-))\\v-v'\ + \hiv)-hiv')\]k,it,x) 

< Clfcl(^,a;)|^;-^;'|^ (7.25) 

for t, v,v' E U,x E M, where ki(t, x) is as in (17.241) . so (17.151) is satisfied. Moreover, 

\fy{t,V,x) - f,o{t',V,x)\ 

= |/i'(t;)|||t-x|'^(^)-^/"log|t-a;| - |t' - xl''^^)-^/" log |t' - x| | 

< k2it,t\x), (7.26) 

for t, t',v E U,x E M, where 

for constants 02,03. Then \\k2{t,t' , < Ci\t — so (17. 16p is satisfied taking rj = h_. 
Strong localisability follows from Theorem 7.3(c). ■ 

To conclude this section we examine stationary moving average processes. These pro- 
vide examples of localisable a-stable processes of a rather different nature being stationary 
processes and not based on existing sssi processes. 

Proposition 7.5 Let < a < 2, let g E Ta o.'nd let M be a symmetric a-stable measure 
on M with control measure C. Define the stationary process Y by 



Y{t) = J 9{t- x)M{dx) {t E R). (7.28) 
Suppose that there exist jointly measurable functions h{t, .) E Ta such that 



lim 



,(r(t + z))-,M_^(^^^^ 



7-7 



dz = ^ (7.29) 



/or all t G M, where 7 + (l/a) > 0. Then F (7 + {1/ a))-localisable at a// n G R with 
local form Y^ = {/ h{t, z)M{dz) : t G R}. 
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Proof. Using stationarity followed by a change of variable z = —x/r and the self-similarity 
of M, 



Y{u + rt) -Y{u) = Y{rt)-Y{0) 

{g{rt — x) — g{—x))M{dx) 



= r^''" j {g{r{t + z)) - g{rz))M{dz) 

where equality is in finite dimensional distributions. Thus 

Yiu + rt)-Yiu) /",, ,^ f f g(r(t + z)) - g(rz) ,\,.,,^ 
^ ^Jyo. -J h{t,z)M{dz)= j [^-^ ^ -h{t,z)jM{dz). 

By [m Proposition 3.5.1] and ^29^, r-^-^/'^(Y{u + rt) - Y{u)) J h{t,z)M{dz) in 
probability and thus in finite dimensional distributions. ■ 

A particular instance of fl7.28p is the reverse Ornstein-Uhlenbeck process, see [TTl 
Section 3.6]. 

Theorem 7.6 (Reverse Ornstein-Uhlenbeck process) Let A > and < a < 2 and let M 
be an a-stable measure on M with control measure C The stationary process defined by 

/CO 
exp(-A(x - t))M{dx) {t e M) 

has a version in D(R) that is 1/a-localisable at all u E with Y^ = L^, where is 
a-stable Levy motion. 

Proof. The process y is a stationary Markov process which has a version in D(M) see 
[TBI Remark 17.3]. It is a moving average process taking g{x) = exp(Ax)l(_oo,o](3;) in 
fl7.28p . It is easily verified using the dominated convergence theorem that g satisfies 
fl7.29p with 7 = and h{t,z) = — l[_jo](-2), so Proposition 7.5 gives the conclusion with 
Fj(t) = -M(H,0])=L,(t). . 



8 Sums over Poisson processes 



In the next section we will set up 'multistable processes', that is a-stable processes where 
a is allowed to vary with t. For this it is convenient to express the random field X{t,v) 
as a sum over a suitable Poisson point process. 

In this section we bring together the basic properties of Poisson sums that we need. 
Let {E,S,m) be a cr-finite measure space. We work throughout with a Poisson point 
process 11 on ii^ x M, with mean measure m x C where C is Lebesgue measure. Thus 11 is 
a random countable subset of ii^ x M such that, writing N[A) for the number of points in 
a measurable A C -E x M, the random variable N[A) has a Poisson distribution of mean 
(m X C){A) with N{Ai), . . . , N{An) independent for disjoint Ai, . . . , A^ C see [H]. 

We define a quasinorm on certain spaces of measurable functions on E. For < a < 
6 < 2 let 

^a,b = J^a,b{E,S,m) = {/ : / is m-measurable with ||/||a,6 < oo} 
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where ^ 

\\f\\a,b = (^j^ |/(x)rm(rfx)) + (^j^ \f{x)\'m{dx)^ . (8.1) 

(Of course || is a norm if 1 < a < 6.) Note that ii a < a' <h' <h then J^a,b C J^a',h' 
and ||/||a',b' < c||/||a,fe where c depends on a, a', 6', 6. Moreover, J^a,a = Ta- 

The following estimate will be useful. Note that expressions such as (18.20 have two 
parts since we need to control the growth of g{x^ y) at both small and large values of y. 

Lemma 8.1 Let g : E x M. ^ be C"^ -measurable and suppose that 

\gix,y)\<hix){\y\-'/^+\y\-'/') (8.2) 

where h G J^a,b for some < a < b < 2. Then there is a constant c depending only on a 
and b such that 




sm^{\eg{x,y))m{dx)dy <c\e'' j \h{x)\''m{dx) + 9^ j \h{x)\^m{dx)j (^ > 0). 

(8.3) 

Proof. We have 

sm^ {^9 g{x,y))m{dx)dy < J J mm \^j9'^\g{x,y)\'^ ,1^ m{dx)dy 

< ci j j mm{9^\h{x)\^\y\-^/'',l}m{dx)dy 

+ci [ [ mm{9^\h{x)\'^\y\-^^\l}m{dx)dy, (8.4) 




where ci is a constant, using (18. 2p and making a simple estimate. But 

min {9^\h{x)\^\y\~^^'', l}m{dx)dy 




< [ [ dy + 9^\h{x)\^ [ \yt^'''dy 

J L"'|i/|<|6»/i(xOI'' J \y\>\dh(x)\°- 

< ca^ / \h{x)\''m{dx) 



m{dx) 



where C2 depends only on a, so along with a similar estimate with b replacing a, (18. 4p 
gives (18. 3p . ■ 

The next proposition gives criteria for the convergence of Poisson sums. We write 
(X, Y) for a random point of -E x R of the Poisson process 11. 

Proposition 8.2 Let g: Ex'R^Wbemx C-measurable with 

\g{x,y)\<h{x){\y\-'/''+\y\-'/'') (8.5) 

where h G J-'a,b- 
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(a) //O < a < 6 < 1 then the series 

5^ ^?(X,Y) (8.6) 
(x,Y)gn 

converges absolutely almost surely. 

(b) Suppose that < a < b < 2 and that g is symmetric in the sense that 

g{x, -y) = -g{x, y) (x, G ^ x M. (8.7) 

Let En be an increasing sequence of m-measurable subsets of E with m{En) < oo for all 
n and U'^^^En = E and write Rn for the rectangle {{x,y) : x G En, \y\ < n} G E x M. 
Then we may define 

Yl ^?(X'Y)= lim Yl (8-8) 

(x,Y)Gn (x,Y)Gnnfl„ 

where the series converges almost surely. 

(c) Provided the symmetry condition (8.7) holds, the characteristic function of Tj, 



taking either definition { \8.6^ or definition { \8.8^ . is given by 

E,e-)^exp(-2//.n^,l.«(..,))™(<^.)*) (« . M). (8.9) 

Proof. If < a < 6 < 1, (18. 5p easily implies that J mm{\g{x,y)\,l}m{dx)dy < oo. 
By Campbell's theorem [12i Section 3.2] the random sum (18.61) is absolutely convergent 
almost surely with characteristic function 

E(e-,^exp(//(e-.)-l)»(..M.) 

If the symmetry condition (18. 7p holds, this reduces to (18.91) . 

In case (b) write S„ = E(x,Y)ennii„ ^(X, Y) = E(x,Y)en^(X, Y)1k„(X, Y), where 1r„ 
is the indicator function of Then by (18. 5p / mm{\g{x,y)lR„{x,y)\, l}m{dx)dy < oo, 
so using Campbell's theorem just as before 



E(e^^^") = exp\^-2j j sm\^9g{x,y))lR„{x,y)m{dx)dy 
— exp I — 2 / / sm^{\9g{x,y))m{dx)dy 



as n — > oo for all 9, by monotone convergence. By (18. 3p there is a number ci > such 
that 



1 > exp -2 / / ?,iYi'a9g{x,y))m{dx)dy > exp(-2ci|e|") > 1 - 2ci|e|'' 



for 16*1 < 1, using that 1 — < x if x > 0. Thus lim^^oo E(e*^^") exists for all 9 and is 
continuous at 6* = 0, so by Levy's continuity theorem [5l Section 10.6], S„ converges in 
distribution to a random variable S with characteristic function (18.91) . 
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We may write 

lim ^?(X'Y) = J] Yl ^(X,Y) 

n— >oo ' ' < ' < ' 

(x,Y)ennfln n=\ (x,Y)6nn(fi„\fl„_i) 

(taking i?o = 0)? which is an infinite sum of independent random variables that converges 
in distribution, so by another theorem of Levy [5], Chapter 12] it also converges almost 
surely. ■ 

Proposition 8.3 Let S = X](x Y)en ^) (1^.61) or fjg.gp where g{x,—y) = 

—g{x,y) and 

\gix,y)\<hix){\y\-'/''+\y\-'/') (8.10) 
for some h G J-'a,b- Then for < p < a, 

E|Sr<c||/^||^,„ (8.11) 

where c depends only on a, b and p. 

Proof. A simple calculation using characteristic functions (see ^ P-47]) gives 



-2/X 



P{|S|>A} < -/ {l-E{exp{i9E)))d9 

2 J-2/X 

--(-//-(^--)M--))^^ 

l-exp(-2c(r/|M.)IM<^.)+^'/lM.)l'™W)))<^^ 



A f'/' 
2 



-2/A 
2/A 

b I li,/„Mf)^ 



< cA y (^ry \h{x)\''m{dx) + e" J \h{x)\''m{dx) j de 

< c.A-/iM.)r,„(..) + c.A-'/iM.)r™(<i.) 

= A-"/i, + A-%, 
say, where Ci depends only on a and 6, using (18.91) and (18. 3p . Then 

/■oo 

E|E|P = p X^-^Pim > X)dX 
Jo 

POO POO 

< p XP-^ mm{l,X''' ha} dX+p / A^-^ min{l, A~%}rfA 

Jo Jo 

< P I XP^^dX + pha / A^-^'-^c/A + p / ' A^-^c/A + p/ib / AP-''-irfA 

io ihy- Jo Jhi'" 

< o,{ht^hf) 

< c\\h\\lb 

where C2, c depend on a, 6 and p. ■ 

We will sometimes need the following variant of Proposition 8.3. 
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Corollary 8.4 Let 0<p<a<ai<bi<b<2 and let /i, /2 G J-'a,b- Let 



^= Yl (/i(X)Y<-i/"i> -/2(X)Y<-^/"2>) , 
(x,Y)en 

where ai < ai, 02 < Then 

m'<c\\h-f2\\% + c\\hrja,-a2\' 
where c depends only on a, ai, b, bi and p. 
Proof. Since 

E = ^ (/i(X) -/2(X))Y<-i/"i> + /2(X) J2 (Y<-i/"i> - Y<-i/"2>) 
(x,Y)en {x,Y)Gn 

= Yl (/i(X)-/2(X))Y<-i/"^> + /2(X) Y («i-«2)Y<-V">a-2log|Y| 
(x,Y)en (x,Y)en 

where a G [ai, 02], using the mean value theorem, the corollary follows from Proposition 
8.3. ■ 

Note that the introduction of ai and 61 in Corollary 8.4 is necessitated by the 'log' 
term to ensure uniformity of the constant c. 



9 Multistable processes 



We now show how our approach may be used to construct multistable processes, that is 
processes where the local stability index varies. The development of this section mirrors 
that of Section 7, but depends heavily on the properties of Poisson sums derived in Section 
8. We seek an analogue of Theorem 7.3 but with the local form an a;(u)-stable process 
with a{u) depending on u. 

We first define a random field analogous to (17.121) . but where the stable random mea- 
sure M is not allied to a particular value of a. Whilst it would be possible to set up a ran- 
dom measure that resembles an a(M)-stable measure close to u, this would be technically 
quite complicated. We therefore favour an alternative approach, using a representation 
by sums over Poisson processes. In particular this permits X{-,v) to be specified using 
the same underlying Poisson process for different v. 

As before {E,S,m) is a cr-finite measure space, and 11 is a Poisson process on x M 
with mean measure m x C In the case of constant a, with M a symmetric a-stable 
random measure on E with control measure m and skewness 0, the stochastic integral 
(17.21) may be expressed ClS db Poisson process sum 

nf) = I f{x)M{dx) = c{a) Y /(X)Y<-i/-> (0 < a < 2), (9.1) 

(X,Y)Gn 

with the sum taken in the sense of (18.61) or (18.81) . and with 

c(a) = (2a^ir(l - a) cos(i7ra))"^^" , (9.2) 
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see [TTl Section 3.12]. (As before a^'^^ = sign(a)|a|'' and C is Lebesgue measure.) 

Particularly relevant in (9.1) is that the stability index a occurs only as an exponent of 
Y, since the underlying Poisson process does not depend on a, so by varying this exponent 
we can vary the stability index. Thus the random field 

X{t,v)= ^ f{t,v,X)Y<~'/"^^> (9.3) 

(X,Y)Gn 

gives rise to a multistable process with varying a, of the form 

Y{t)^X{t,t)= /(t,t,X)Y<-i/-W>. (9.4) 

(X,Y)Gn 

We first consider continuity and boundedness of the processes. 

Proposition 9.1 Let U be a closed interval. Let X be the random field defined by 

X{t,v)= J2 fit,v,X)Y<~'/"^^> {t,veU) (9.5) 
(x,Y)en 

where f{t,v,-) G J-'a,b (n^^ jointly measurable and a : U ^ {a,b) is continuous. 
(a) Suppose < a < a{v) < b < 1 for v E U. If 

sup \fit,v,x)\<k{x), (9.6) 
t,veu 

where k G J-'a,b! then {X{t,v) : t, f G U} has a bounded version. 

If in addition {f{t,v,x) : x G E} is an equiuniformly continuous family for t,v G U, 
then X has a continuous version. 

(6) Suppose that 1 < a < a{v) < b < 2 for v E U and 1/a < rj < 1. Suppose that 

\a{v)-a{v)\<ki\v-v'\'^ {v,v'eU), (9.7) 

that 

sup ||/(t,t;,-)||a,b < oo, (9.8) 

t,veU 

and 

\\f{t, V, ■) - fit', v', ■) ||a,6 <k2i\v- vT + \t- tr) (t, t', V, v' G U) . (9.9) 

Then Y = {X{t,t) : t G U} has a continuous version satisfying an a.s. (3-Hdlder condition 
for alio < (3 < {ija - l)/a. 

Proof (a) From ( 1931) . for all t,v eU, 

\X{t,v)\ < J2 l/(t,^',X)||Y|-V-W 
(x,Y)gn 



< Yl sup |/(t,t;,X)|(|Y|-i/" + |YrV^)^Z 



(x,Y)en*'^^^ 

where Z is an a.s. finite random variable, by Proposition 8.2(a). Thus {X{t, v) : t,v E U} 
is a.s. bounded. 
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Assuming also the equicontinuity condition, given e > we may choose r > 1 such 
that ||/i;(x)l{|2^|>r}(x)||a,6 < e, where 1 is the indicator function. By equiuniform continuity 
we may find 6 > such that for all x G M and \(t,v) — (t',v')\ < S we have \f(t,v,x) — 
f{t'^v',x)\ < r~^^"'e, and \a{v) — a{v')\ < e. Then if \{t,v) — {t',v')\ < 5, making several 
estimates in the obvious way, 

\X{t,v) - X{t',v')\ < Yl |/(t,^^,X)Y<-i/"(^)> -/(t',i;',X)Y<-i/"(''')>| 

(x,Y)en 

< \fit,v,X) - f{t',v',X)\\Y\-'/^^^^ +2^ sup \f{t,v,X)\\Y\-'/^^^^ 

|X|<r |X|>r*'''^^ 

+ Y |/(t',?^',X)||Y<-i/"(^)> -Y<-i/"(^')>| 
(x,Y)en 

< E r-^/'^el||.|<.}(X)|Yri/"(^)+2 MX)l{|.|>n(X)|Yr^/'^('') 
(x,Y)en (x,Y)en 

+ Y l/(t',^',X)|i^|Y|-V"|log| 
(x,Y)en 



aiv] — a[v 



< ( J2 r-'/''el{|.|<,}(X) + 2 J2 MX)l{|.|>n(X) 
(x,Y)en (x,Y)en 

+ E ci|A:(X)|le)(|Yri/'^ + |YrV'')^Z, 



(x,Y)en 

where is a random variable, and we have used the mean value theorem in the third 
term of the sum with a G [a{v),a{v')]. Fix < J9 < a. By (IS.lip there is a constant c 
independent of e such that 

E\Z,\P< ceP. 

The proof is completed just as in the proof of Proposition 7.1(a). 
(b) We estimate 

X{t, v) - X{t', v') = (X(t, v) - X{t, v')) + v') - X{t', v')) (t, t', V, v' G U) (9.10) 
by considering its two parts in turn. Firstly 

X(t, v) - X(t, v') = J2 (/(^' ^' X)Y<-V"('')> - /(t, v', X)Y<"V"K)>) 
(x,Y)en 

Thus Corollary 8.4 gives, for < p < a, 

E\X{t,v)-X{t,vT < ci||/(t,t;,-)-/(t,t^',-)lll. + ci||/(t,t;',-)lllJ«(t;)-a(t;')r 

< C2\v-v'\''P (9.11) 

by dH]), dH]) and 

For the second term of fl9.10p 

Xit,v)-Xit',v)= Y (/(t,^,X)-/(t',^,X))Y<-V"(-)>. 

(X,Y)6n 
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Then 

\{f{t,v,x) - fit\v,x))y<-'/'^^^>\ < \fit,v,x) - fit',v,x)\ {\y\-'/^ + \y\-'/') 
so, for < p < a, Proposition 8.3 and ( 19. 9p give 

E\Xit,v)-Xit',v)\^ < cs\\fit,vr)-fit',vr)\\l, 

Combining with flQ.lip we estimate flQ.lOp to get, for t,t',v,v' G U, 
E\X{t,v) - X{t',v')f < cr,{\v-v'\''^ + \t-t'\'^P). 

Speciahsing, 

E|y(t) -r(t')r = E|x(t,t) -x(t',t')r < 2c5\t-t'\''p 

for t, t' G U. 

Since rj > 1/a we may choose < p < a such that rjp > 1. Kolmogorov's Theorem 
5.2 gives that {Y{t) : t E U} has a continuous version that is a.s. /?-Holder for all 
< (3 < {rip — l)/p for all p < a. m 

We come to the main result on the localisability of processes with varying stability 
index. 

Theorem 9.2 Let U be a closed interval with u an interior point and let < a < b < 2. 
Let X be the random field defined by 

X{t,v)= f{t,v,X)Y<-'/^^^> {t,veU) (9.12) 

(X,Y)Gn 

where f(t,v, ■) G J^a,b o'^e jointly measurable and a : U ^ {a,b). 

(a) Suppose X{-,u) is h-localisable at u for h > 0. Suppose that sup^^jj \\f{t,u,-)\\a,b < 
oo, and that for some rj > h 

\a{v) - a{u)\ < h\v - ul'^ (veU), (9.13) 

and 

||/(t, V, ■) - fit, u, ■) |U„fe <h\v- it, veU). (9.14) 

Then Y = {Xit,t) : t G U} is h-localisable at u with local form y^'(-) = X'^i-,u). 

(6) Suppose that < «(«) < 1 and that Xi-,u) is strongly h-localisable in C(]R) (resp. 
D{M) ) at u. Suppose that for some rj > h 

\aiv) - aiu)\ < ki\v - ul'^ ivEU), (9.15) 

and 

\fit,u,x)\<k2ix) iteU,xeE), (9.16) 

and 

\fit,v,x)- fit,u,x)\<hix)\v~u\'' it,veU,xeE), (9.17) 

where k2i-),k^i-) G J^a,b- If Y = {Xit,t) : t E U} has a version in C(f/) (resp. DiU)) 
then Y is strongly h-localisable in C(]R) (resp. -D(R)^ at u with Y^i-) = X^(-,m). 
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(c) Suppose that 1 < a{u) < 2 and that is strongly h-localisahle in C(M) (resp. 

D(M)j at u. Let rj satisfy l/a{u) < rj < 1. Suppose that a is continuously diff'erentiable 
on U with 

\a\v) - a'{v')\<ki\v -v'\^ {v,v'eU). (9.18) 

Suppose that the partial derivatives fy{t,v,-) G J^a,b for all t,v & U, and the following 
estimates hold: 

sup \\f{t,U, ■)\\a,b < OO, sup \\fy{t,V,-)\\a,b < OO, (9.19) 

t(^u t,veu 

||/(t, V, ■) - fit', V, ■) ||a,fe <k2\t- t'r (t, t', veU), (9.20) 

\f,it,v,x)- f,{t,v',x)\<ks{t,x)\v-v'\'^ {t,v,v' eU,xeE), (9.21) 

and 

\f,{t,v,x)- f,it',v,x)\<ki{t,t',x) it,t',veU,xeE), (9.22) 

where sup^^^ Wk^iJ:, ■)\\a^f, < oo, and ||fc4(t, t', OlU.fe < k\t — t'\^ for all t,t' G U. Then 
Y = {X(t,t) : t E U} is strongly h-localisahle in C(R) at u with FJ(-) = X^(-,u). 

Proof, (a) We have 

X(t, v) - X{t, u)= J2 (/(^' X)Y<-i/"('')> - /(t, u, X)Y<-i/"(")>) . (9.23) 
(x,Y)en 

With < p < a, Corollary 8.4 gives that there are constants Ci, C2 such that 

E\X{t,v)-X{t,u)\^ < c4f{t,v,-)-f{t,ur)\\% + cMv)-a{u)nf{t,u,-)\\l, 

< C2\V-U\''P, 

for all t,v E U, by fl9.13p and (19.141) . Part (a) now follows from Theorem 4.1. 

(b) We may assume that a < a{v) < 6 < 1 for t> G f/, if necessary using the continuity 
of a to replace [/ by a subinterval to decrease the value of b. Splitting and estimating 
fl9.23p . using the mean value theorem as in the proof of Corollary 8.4, we get, with 
a' G [a{v),a{u)] (where a' depends on v), 

\X{t,v)-X{t,u)\ < \fit,v,X) - f{t,u,X)\\Y\-'/"^^^ 

(X,Y)Gn 

+ Yl \f{t,u,X)\\a{v)-a{u)\\Y\-'/^'a'-'\\og\Y\\ 
(x,Y)en 

< It; -Ml" J2 l'^3(X)|(|Y|-i/'^ + |Y|-i/') 

(X,Y)Gn 

+ci|i;-nr Yl l^2(X)|(|Y|-i/'^ + |Y|-^/^) 
(x,Y)en 

for all t,v eU, using fl9J[5|) - fl9Tri) . By Proposition 8.2(a) |fc2(X)|(|Y|-i/" + \Y\~^/^) and 
|^3(X)|(|Y|~-'^/'* + lYI^-*^/^) are a.s. finite random variables, so (15.21) holds and Theorem 5.1 
implies that Y is strongly localisable at u. 

(c) The conditions of Proposition 9.1(b) are easily checked, so Y has a continuous 
version. We may assume that 1 < a < a{v) < b for v E U and that 1/a < 77, using 
continuity of a to replace ?7 by a subinterval and to increase the value of a if necesssary. 
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Define 



, , X(t,v) - X(t,u) , , , 



V — u 



again we use Kolmogorov's criterion to show that {Z{v,v) : v G U} is almost surely 
bounded to get (15.21) . We write 



Z{t,v) = Zi{t,v) + Z2{t,v) {t,veU) 



(9.24) 



where 



Z,{t,v)= J2 9{t,v,X)Y<-'/^^-> with g{t,v,x)= ^^^''''''] fit,u,x) ^^^5) 



(x,Y)en 



V — u 



and 



{X,Y)Gn 



Y<-l/a(f)> _ Y<-l/a{u)> 



V — u 



(9.26) 



For p < a we estimate E|Z(t, v) — Z{t', v')\^ by breaking it into four parts. 

(i) Applying Lemma 7.2 to (I9.2ip gives \g(t,v,x) — g{t,v',x)\ < 2^k3{t,x)\v — v'l"^. 
Thus Corollary 8.4 on (I9.25P and then fl9.19p with the mean value theorem gives 

E\Z^{t,v)-Z,it,v')r < ci||^(t,t;,0-<7(t,t;',-)lll. + ci|k(t,t;',-)li:j«(^)-«(^')r 

(9.27) 



< c,\Mt,-)\\l,\v-vT + C2snp ||/.(t,., 

t,veu 



(ii) Using the mean value theorem and ( 19.22^ 

1 



\g{t,v,x) -g{t',v,x)\ 



\{f{t,v,x) - f{t',v,x)) - {f{t,u,x) - f{t',u,x))\ 



\v — u\ 

= \fy{t,Vi,x) - f^{t',vi,x)\ 

< hit,t',x) 



where Vi G {v,u) depends on t,t',v and x. Proposition 8.3 with (19.250 now gives 

E\Z,{t,v) - Z,{t\v)\P < c,\\h{t,t'r)\\l, < c,\t-t'r. (9.28) 
(iii) Turning to (I9.26p . a simple estimate using (19.181) gives that 



^y<-l/aiv)> 

dv 



< Cq\v — v' 



By Lemma 7.2 



^<-l/a(D)> _ y<-l/a{u)> y<-l/a(v')> _ y<-l/a(u)> 



V — U 



+ \y\-n- 



<2^c,\v-vr {\y\-"''+\y\-"'). 



Thus Proposition 8.3 applied to (I9.26P gives 



E|Z2(t, v) - Z2(t, v') r < cr||/(t, ■) ||^,> - v'r <cs\v- v' 



VP 



(9.29) 
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(iv) Finally, a further mean value estimate gives 

<~l/cx(v)> _ <-1/q(«)> 



V — u 

< c,\{f{t,v,x) - f{t',v,x))\ {\y\'/'^ + \y\'/') . 

By Proposition 8.3 and fICTD 

E\Z2it,v) - Z,{t\v)Y < cr,\\f{t,vr)-f{t',v,-)\\l, < c,,\t-t'r. (9.30) 

Taking fICTD with dOTj) . (lOSD . (19:291) and (ICTD . we conclude that for some Ci2 
independent of t, t', t;, f ' G t/, 

E|Z(t,t;) -Z(t',t;')r < cisd^^ - ^^T^ + |t - tT^) (9.31) 
if < p < a. Specialising, 

E\Z{v,v)-Z{v\vX < 2ci2K'-t/r^ (9.32) 

for V, v' e U . 

Since r/ > 1/a we may choose < p < a such that rjp > 1. Applying Kolomogorov's 

Theorem 5.2 to {Z(^;,^;) : v E U} we conclude that Z(v,v) = ^ ^ — has a 

f — u 

version that is a.s. bounded. Thus (15.21) holds and strong localisability follows from 
Theorem 5.1. ■ 

We now show how Theorem 9.2 may be used to construct some specific multistable 
processes. In these examples we take {E,S, m) to be Lebesgue measure on M so that 11 
is the Poisson process on with mean measure C^. We first construct a multistable 
analogue of the linear multifractional motion of Theorem 7.4. 

Theorem 9.3 (Linear multistable multifractional motion). Let a : M — be and 
a -.R-^ (0,2) and h -.R^ (0, 1) be C^. Define 

Y{t) = a{t)c{a{t)) Y<-i/"(*)> {\t - X\^(i)-y^(t) _ |x|Mt)-i/"W) ^ (9.33) 

(X,Y)Gn 

(a) The process Y is h{u)-localisable at all u G M, with = a{u)Lo,(u),h{u), where L^^h 
is linear stable motion. 

(6) If u is such that h{u) > l/a{u) then Y is strongly h{u)-localisable in C(M) at u, 
with Y^ = a{u)La{u)M'^)- 

Proof. By the amplitude result. Proposition 3.2, the term a(t)c{a{t)) in (19.331) does not 
affect localisability, so it is enough to prove the result with a(t)c{a{t)) = 1. Define a 
random field by 

(x,Y)en 

= Yl /(i,^,X)Y<-V"(^)>, 
(x,Y)en 
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where 

f{t,v,x) = {\t - xi'^^^)-^/"^^) - \x\''^^yy'^i^^) . 

Then 

fv{t,V,x) 

= {\t - a; I V^C'') log \t-x\- |a;|^(^)-V"{") log {h'{v) + a'{v)/a{vf) . 

Given m e M we may use continuity of h and a to choose U to be a small enough closed 
interval with u an interior point, and numbers a, b, h-, h+, such that < a < a{v) < b < 2 
and < h- < h{v) < h+ < 1 for all v eU, and such that ^ - | < /i_ < /i+ < 1 - - |). 
A similar argument to that of Theorem 7.4 gives that 

\f(t,v,x)\,\Ut,v,x)\<ki(t,x) (t,veU,xeR) (9.34) 

and 

\f{t,v,x)-f{t,v',x)\, \U{t,v,x)-U{t,v',x)\ < h{t,x)\v-v'\ {t,v,v' eU,xeR) (9.35) 
where 

I. u N / cimax{l,|t-x|'^— + V«} (|x| < 1 + 2 max^et; |t|) , . 

f^,[t,x) - I ^1^1 > i + 2maxiet/|t|) ^ ^ 

for appropriately chosen constants Ci and C2. By virtue of the conditions on a,b, h_, 
it follows that sup^gt; ||/ci(t, •)||a,(, < oo. Since X{-,u) = c{a{u))~^La{u),h{u){')j Theorem 
9.2(a) gives /i(M)-localisability of Y with Y^{-) = X;(-,m) = c{a{u))-\L^(^u)Mn))'u{-) = 
c(a(u))~^La(u)Mu)(')- 

For part (b), we choose U and the numbers a, b, to satisfy the conditions stip- 

ulated in the proof of (a) but also to satisfy > l/a + (l/a — 1/6) > 0, so in particular 
h{v) — l/a{v') > for all v,v' G U. Again as in Theorem 7.4, 

\f{t,v,x) - f{t',v,x)\, \fy{t,v,x) - fy{t',v,x)\ < k2{t,t',x) 

for t,v E U,x eR, where 

l^2[t,t,X)-<^ C4|x-i(i-i')|'^+-'/'-'|^-^'| i\x-l{t-t')\>\t-t'\) 

for constants C3, C4. Then ||A;2(^, ■)\\a,b < c^lt — The conditions of Theorem 9.2(c) 

are satisfied with rj — 1/a > l/a{u), so strong localisability follows. ■ 

Note that the differentiabihty conditions in Theorem 9.3 could be weakened slightly 
to Holder conditions for which Theorem 9.2 would still be applicable. 

Recall that an a-stable Levy motion, < a < 2, is a process of D(R) with station- 
ary independent increments which have a strictly a-stable distribution. Taking M as a 
symmetric a-stable random measure on R, the a-stable Levy motion may be represented 
as 

L„(t) = M[0,t] = I l[o,t](a;)M(cia;) = c(a) l[o,^](X)Y<-V"> {t e M), (9.38) 

{X,Y)Gn 

where 11 is the Poisson process on with as mean measure, l^^t] is the indicator 

function and c(a) = (2a^^r(l — a) cos(|7ra)) Then La is 1/a-sssi and is strongly 
1/a-locahsable in D{ 
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Theorem 9.4 (Multistahle Levy motion). Let a : ^ (0,2) and a : M ^ M+ &e 

continuously differentiahle, and define 

Y{t) = a{t)c{a{t)) l[o,t](X)Y<-i/"W> (t G R). (9.39) 

(x,Y)gn 

(a) If I < a{u) < 2 then Y is 1/ a{u)-localisahle at u, with Y^ = a{u)La(u)- 

(h) If < a{u) < 1 and a'{u) ^ then Y is 1-localisable at u with {Y^^it) : t G M} = 
{tW : t G M}, where W is the random variable 

W = a{u)c{a{u)) J] „j(X)Y<-V"W> log |Y| | + |-(a(«)c(«(t.)))y 

(x,Y)en \ \ > J 

Proof, (a) By Proposition 3.2 the term a{t)c{a{t)) does not affect locahsabihty. Define a 
random field by 

X{t,v)= Yl l[o,t](X)Y<-^/'^(^)> {t,veR). 

(X,Y)Gn 

Taking f{t,v,x) = l[o,t](a;) the conditions of Theorem 9.2(a) are satisfied with h = 
l/a{u) < 1 = rj, so the result follows from the localisability of L^. 

(b) In the case where a{t)c{a(t)) = 1 



Y{u + ri) Y{u) = i ^ j^^^^(^x^)(^Y<-i/a(«+rt)> _ Y<-i/aM>^) 
^ ^ (x,Y)en 



/<-l/a{u+rt)> 

r 

(x,Y)en 



Letting r — >■ the second term vanishes if l/a{u) > 1 and the first term converges to W 
in finite dimensional distributions. The general case is similar. ■ 

Note that Theorem 9.4(b) illustrates a general phenomenon that occurs when the 
process {X{t,u) : t G M} is /z(M)-localisable at u where h{u) > 1. The process {Y^{t) : 
t G M} will typically be 1-localisable at u, with the dominant component of Y^{t) derived 
from {X{u + rt,u + rt) — X{u + rt, u))/r rather than from X'^{t, u). 

As explained in [17, Section 7.6], there are two ways to extend the linear fractional 
stable motion to the case H = 1/a. Apart from the Levy motion considered above, one 
may define the following process, called log-fractional stable motion: 

/oo 
{\og{\t-x\)-\og{\x\))M{dx) (tGR) (9.40) 
-oo 

where, as usual, M is an a-stable random measure. This process is well-defined only 
for a G (1,2] (the integrand does not belong to JF^, for a < 1). It is 1/a-self-similar 
with stationary increments. Unlike the Levy motion, however, its increments are not 
independent. Another difference is that log-fractional stable motion does not have a 
version in -D(M), so we cannot speak of strong localisability. 
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Theorem 9.5 (Log-fractional multistable motion) Let a : M — > (1,2) and a be continu- 
ously differentiahle, and define 

Y{t) = a{t) (log|t-X| -log|X|)Y<-i/"(*> (teM). (9.41) 

(x,Y)en 

Then Y is 1 / a{u)-localisable at all ti e M, with Y^ = a{u)Aa(u)- 

Proof. The proof is similar to that of Theorem 9.4, by considering the field 

X{t,v)^ (log|t-X| -log|X|)Y<-^/"(''> (^,^;eK), (9.42) 

(X,Y)Gn 

with Theorem 9.2(a) is applied to f{t, v, x) — log \t — x\— log \x\. ■ 
For a final example we give a multistable version of Theorem 7.6 

Theorem 9.6 (Multistable reverse Ornstein-Uhlenbeeck process) Let A > and a : M — >■ 
(1,2) be continuously dijferentiable. Let 

Y{t)^ Y exp(-A(X - f))Y<-V«(*)> (feR). 
(x,Y)en,x>t 

Then Y is 1 / a{u)-localisable at all u & W, with Y^ — c{a{u))~^La{u)) where L^ is Levy 
a-stable motion. 

Proof. Let 

X{t,v)= Y exp(-A(X - t))Y<-^/"(^> (t,weM). (9.43) 

(X,Y)Gn,X>t 

Then for each v, X(-, v) is everywhere l/a(t')-localisable with X^(-, v) = c(a(t'))~^La(^)(-) 
by Theorem 7.6. Applying Theorem 9.2(a) with f{t,v,x) — l[t,oo) exp(— A(x — t)) gives 
the conclusion. ■ 

10 Further work 

There are a great many possible variants and extensions of this work. Localisable pro- 
cesses of many other forms may be constructed. For example multistable processes with 
skewness and the class of stationary localisable processes deserve investigation. There 
may be advantages in seeking other representations of multistable processes such as by 
sums involving arrival times of a Poisson process or as stochastic integrals with respect to 
suitably constructed random measures. Our conditions for localisability could certainly 
be weakened and further techniques for establishing localisability and in particular strong 
localisability developed. It would be interesting to study long range dependence of Mul- 
tistable processes. Effective techniques for simulation and inference on parameters for 
these processes are also needed. We will be addressing some of these matters in a sequel 
to this paper. 
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